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Abstract 

Recently, observational searches for gravitational wave background (GWB) have been developed 
and given direct and indirect constraints on the energy density of GWB in a broad range of 
frequencies. These constraints have already rejected some theoretical models of large GWB spectra. 
However, at 100 MHz, there is no strict upper limit from direct observation, though the indirect limit 
by ^He abundance due to big-bang nucleosynthesis exists. In this paper, we propose an experiment 
with laser interferometers searching GWB at 100 MHz. We considered three detector designs and 
evaluated the GW response functions of a single detector. As a result, we found that, at 100 MHz, 
the most sensitive detector is the design, a so-called synchronous recycling interferometer, which 
has better sensitivity than an ordinary Fabry-Perot Michelson interferometer by a factor of 3.3 
at 100 MHz. When we select the arm length of 0.75 m and realistic optical parameters, the best 
sensitivity achievable is /i 7.8 x 10"^^ Hz"^/^ at 100 MHz with bandwidth ~ 2 kHz. 

PACS numbers: 04.80.Nn, 95.55.Ym 
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I. INTRODUCTION 



There are many theoretical predictions of gravitational wave background (GWB) in a 
broad range of frequencies, 10^^^ — 10^° Hz. Some models in cosmology and particle physics 
predict that there are large stochastic GWB at ultra high frequency ~ 100 MHz. Testing 
these models with GW detectors for high frequencies is very important. In the quintessen- 
tial inflation model [1], the blue spectrum of GWB is produced during the kinetic energy- 



dominated era after the inflationary expansion of the universe 
to subsequent reheating processes 



and is very sensitive 



^. In other inflation models, during the first stage of the 
process of reheating, called preheating, GWB at high frequencies is created due to large den- 

8j. Pre- big-bang scenarios in string cosmology can also generate 



sity inhomogeneities 

m 

high-frequency backgrounds 
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Not only cosmological sources but also compact ob- 
jects can create GWB at ~ 100 MHz. Primordial black holes produced in the early universe, 
which have much smaller masses than the sun, emit GW via binary evolution and coales- 
cence 



12 



131] and evaporation [14|. Recent predictions of GW emission from black strings 



in the Randall-Sundrum model generate spectral features characteristic of the curvature of 



extra dimensions at high frequencies 15|, [16 1 . 



Upper limits on GWB in wide-frequency ranges have been obtained from various obser- 
vations; cosmic microwave radiation at 10^^^ — 10^^^ Hz 17|, pulsar timing at 10^^ — 10"^ Hz 
[isl . Doppler tracking of the Cassini spacecraft at 10~^ — 10~^Hz [iqI], direct observation 
by LIGO at 10 — lO'^Hz 20| . ^He abundance due to big-bang nucleosynthesis at greater 
frequencies than 10~^°Hz |21[]. Nevertheless, as far as we know, no direct experiment has 



22|. 



been done above 10^ Hz except for the experiment by A. M. Cruise and R. M. J. Ingley 
They have used electromagnetic waveguides and obtained an upper limit on the amplitude 
of GW backgrounds, h < 10~^^ corresponding to hlQ^Qgy, < 10^'^ at 100 MHz, where /iioo is 
the Hubble constant normalized with lOOkmsec"^ Mpc~^ and Qg„ is the energy density of 
GWB per logarithmic frequency bin normalized by the critical energy density of the universe, 
that is fig„(/) = {dpgfj / d\n. f) / pc- This constraint is much weaker than the constraints at 
other frequencies. Therefore, a much tighter bound above 10^ Hz is needed to test various 
theoretical models. 

In this paper, we propose a method of direct detection of GW at 100 MHz with laser inter- 
ferometers. At high frequencies, the GW wavelength is comparable to the size of a detector. 
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which is the order of a few meters. The usual approximation that the GW wavelength is 
much larger than the detector size is not valid. In this case, the phase of GW changes during 
the one-way trip of light between mirrors. Therefore, we have to use a detector design that is 
able to integrate GW signals efficiently. Note that our investigation is general and applicable 
to other than 100 MHz. The contents of the paper are as follows: In Sec. II, we will consider 
three detector designs and give GW response functions. Detailed calculations are given in 
the Appendices. In Sec. Ill, the GW responses numerically calculated are compared. Sec. 
IV is devoted to discussions and conclusions. 



II. DETECTOR DESIGNS AND THE RESPONSE FUNCTIONS 

In this section, we will consider three detector designs, (i) synchronous recycling interfer- 
ometer (SRI), (ii) Fabry- Perot Michelson interferometer (FPMI) and (iii) L-shaped cavity 
Michelson interferometer (LMI), and derive the response functions for GWs. 



A. Detector designs 

SRI (Fig{T]) was first proposed by R. W. P. Drever in |23| and detailed calculations 



have been done in |24j . |25| |. Laser light is split at a beam splitter and sent into an SR 
cavity through a recycling mirror, which is mirror A located at Xi in Fig{T] The beams 
circulating clockwise and counterclockwise in the cavity experience gravitational waves and 
mirror displacements, leave the cavity, and are recombined at the beam splitter. Then, the 
differential signal is detected at a photodetector. The advantage of SRI is that GW signals at 
certain frequencies are accumulated and amplified because the light beams experience GWs 
with the same phases during round trips in the folded cavity. Consider GW propagating 
normally to the detector plane with an optimal polarization. In this case, the GW signal is 
amplified at the frequencies / = {2n — 1) x c/4L, n = 1,2, ■ ■ ■ , where c is the speed of light 
and L is the arm length. More precisely, the arm length is the distance between Xi and X2 
(or X3) in FiglH On the other hand, the disadvantage of SRI is less sensitivity for GWs at 
low frequencies, / < c/4L, because the GW signal is integrated in the cavity and canceled 
out as the frequency is low. 

The competitive design of detectors with SRI is an ordinary FPMI (Figl2]). FPMI is 
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the most popular design for current ground-based interferometers [26|, Wl, l28|, l29| since it 
has good sensitivity at low frequencies due to the amplification of GW signals with Fabry- 
Perot cavities. However, to amplify the GW signals at high frequencies, one needs to use 
resonance due to the cavity. FPMI has the resonance of GW signals at the frequencies, 
f = n X c/2L, n = 1,2, ■ ■ ■ when the GW response is averaged over the entire sky. To take 
advantage of the resonant response to GW at lOOMHz, the arm length of FPMI should be 
1.5 m. With this detector, one can achieve good sensitivity with narrow bandwidth as well 
as SRI. 

Another possible design of detectors is LMI (Figl3]), whose optical configuration is the 
same as the L-shaped FPMI. However, the GW response resembles SRI rather than FPMI. 
Thus, this design can be regarded as being intermediate between SRI and FPMI. 

To compare the detectors for GWB at 100 MHz, it is necessary to derive the GW response 
functions for GW propagating in arbitrary directions and to compare those averaged on 
the celestial sphere. For GW propagating normal to the detector plane with an optimal 
polarization, it is trivial that SRI and LMI with arm length 0.75 m have maximal sensitivity, 
while FPMI with 1.5 m is not sensitive at all, at 100 MHz. However, FPMI has nonzero 
sensitivity for GW not propagating orthogonally to the arms of the detector since the GW 
response of light going and coming differs. Furthermore, the geometries of detectors also 
affect the GW responses. Therefore, it is nontrivial which is the most sensitive detector. 
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FIG. 1: Synchronous recycling interferometer (SRI). 
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FIG. 2: Fabry- Perot Michelson interferometer (FPMI). 
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FIG. 3: L-shaped cavity Michelson interferometer (LMI). 



B. GW response functions 



The general expression of the phase shift of light induced by GW propagating in an 
arbitrary direction has been derived in many references, for example, 30|, l3l||- We will use 
the expression in [s^ as our starting point. When light travels between two test masses 
located at Xj and Xj, the phase shift of light created by GW in TT (transverse-traceless) 
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gauge is expressed as 

^ r-Lij/c 

5(t)ij{t) = {iiij (g) Yiij) : - dt'y2 ePhp[t - {Lij + ■ X,)/c + (1 - ■ nij)t'], (1) 

Jo p 

where t is the reception time of hght at Xj, uj is the angular frequency of hght, hp is the 
amphtude of GW with plus or cross polarization and is a unit vector in the direction of 
GW propagation. The arm length and unit vector in the direction of the arm are defined 
as Lij = |Xj — Xj| and n^j = (Xj — 'Ki)/Lij, respectively. The symbol : means contraction 
between a tensor and vectors. The tensors e^, p = +, x are polarization tensors of GW and 
are defined as 

e+ = e^^e^-By^ey , (2) 

Gx = Gx Gy -\- Gy iS) Gx 1 (3) 

where e^, and Gy are the unit vectors, which form the orthogonal coordinate with e^. We 
assume that there is no displacement noise, for example, thermal noise, seismic noise, radi- 
ation pressure noise, etc., at 100 MHz. (In fact, rough estimates show that this assumption 
is at least valid as long as the detector sensitivity is above h ~ 10^^'^ Hz^^^^. In the exper- 
iment that reaches better sensitivity, thermal noise of mirrors might limit the sensitivity, 
though other noises are far below. In that case, one should calculate the noise spectrum 
accurately with numerical simulation including the effect of the mirror's complicated re- 
sponse at 100 MHz, which is beyond the scope of this paper.) Under the assumption of no 
displacement, the positions of mirrors in the absence of GW are not perturbed and are just 
given by Xj and Xj. The Fourier transform of Eq. (fl]) is given by 

^-iQ{ez-Xj)/c _ ^-iQ{Lij+ez-Xi)/c 

6c^^,{n) = n., ® n., : - G^hp .^(^ _ . , (4) 

where hp is the Fourier component of the GW amplitude, and Q is the angular frequency of 
GW and is related to the GW frequency with Q = 27r/. 

In general, the response function of a detector is represented by the round-trip signal in 
cavities multiplied by an amplification factor in cavities. We denote the phase shift of the 
round-trip signal by 6(f){Q) and the amplification factor by a{Q). Then, the total output from 
the detector is written as ^^(fi) = Detailed calculations are described 

in the Appendices. We show here only the results. Note that we change the notation of 
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the unit vectors directed in arms and the reflectivity of mirrors in order to simplify the 
expression and make it easy to compare. The response functions of each detector are 

S'^aiim = ain, Rp, RE)S\ui^) , (5) 



^ ' ' ' {Re - Re){1 - RfRe e-^'^^) ' 

and 

• SRI 

Replacing ni2 — > u, 1113 — > v and {Ra, Rc) (Re, Re) in Eq. flAlSI) . the response 
function is 

p 

■ I — J i sinQr + (e, ■ v)(e-*^"'=^-^ - cos fir)} 



ReT, 



1- -v 



U (X> u 

1 - fe, ■ u 



{i sin fir + (e^ ■ u)(e-^''™^-" - cos fir)} ] , (7) 



FPMI (doubled the arm length, r 2r) 

Replacing ni ^ u, n2 — > v in Eq. ( IBlOp . the response function is 



in 



^ { t sin2fir + (e, ■ u)(e-''^™-" - cos2fir) } 



1 — (e^ ■ 

- { i sin2fir + (e, ■ v)(e-"*''™^-^ - cos2fir) } 



1 - (e. ■ v) 
LMI 

Replacing rii — > u, n2 ^ v in Eq. ( ]C8I) . the response function is 



-'all 



(fi) = |^^e^/,pe-^^(^^+''-^-/^) 

p 



^ { i sin fir (e~*Pi + e'^^^ - 2 cos fir) 



l-(e.-v)2 

+ (e, ■ ^){e-'^Pi+P2) + cos2fir - (e'^^^ + e-'""^) cos fir) } 

U (g) u 



{ i sin fir (e-*^^ + e'^^^ - 2 cos fir) 



1 - (e. • u)2 

+ (e, • u)(e~^(Pi+P2) + cos2fir - (e'^^^ + e-'^^) cos fir) } ] . 

(9) 



The phases pi and p2 are defined by pi = Qr (e^ ■ u) and p2 = fir (e^ • v) , r is defined by 
T = L/c and Rp and Re are the amphtude refiectivities of front and end mirrors of cavities, 
respectively. Tp is the amphtude transmissivity of a front mirror of cavities. Note that, in 
the case of SRI, the front and end mirrors correspond to a recychng mirror and three other 
mirrors. Here we doubled the arm length of FPMI so that the first resonant frequency of 
GW signal coincides with that of SRI and LMI. 

III. DETECTOR COMPARISON 

As mentioned in the previous section, in general, the GW response function has the form 
S^aii = Oi We will consider a, which is the common factor for all detectors, and 6(f), 

which depends on the geometry of each detector, separately. 

From Eq. ([6]), the magnitude of the optical amplification factor in the cavities is deter- 
mined only by the (composite) reflectivities of the front and end mirrors. The frequencies 
of the peaks depend on the arm length of detectors, a is plotted in FiglH In the figure, 
we selected L = 0.75 m so that the first resonant peak is located at 100 MHz. At higher 
frequencies, there are many resonant peaks. At lower frequencies, optical amplification is 
stronger as the frequency is lower, since the wavelength of light is larger than the arm length 
of a detector. Here we selected the amplitude refiectivities Rp = 0.99 and Rp = 1 for an 
illustrative purpose. However, in a real experiment, one should select the refiectivities of 
the front mirrors much higher in order to achieve better sensitivity, though the bandwidth 
becomes narrower. 

To evaluate the round-trip phase shift due to GWs Scp^n, we introduce coordinates here. 
The detectors are located on the X-Y plane. Two unit vectors u and v are written as 
u = (0, 1, 0) and v = (1, 0, 0), respectively. We denote the direction of GW propagation by 
the unit vector e^, and the two unit vectors normal to and orthogonal to each other by 

and ey. These are written as 

e^; = ( cos 6 cos 0, cos 6 sin 0, — sin 6 ) , 

ey = { — sin (j), cos 0, ) , 

62 = ( sin 6' cos 0, sin 9 sin 0, cos 9 ) . 

63; and By define the GW polarization tensor in Eqs. ([2]) and ([3]). Here we normalize and 
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FIG. 4: The optical amplification factor a{^l). Parameters are selected L = 0.75m, Rp = 0.99 and 
Re = 1. 



redefine the GW response function as a dimensionless response function, namely, 

r(fi,0,^,^)^i^, (10) 

[UjhT) 

where we assumed that GW has the form ^GpW = h{e^ cos2iIj + ex sin2ilj). ip is the 
polarization angle of GW. Integrating this function by (p, 6 and ip on the celestial sphere 
and averaging lead to 

rLsm ^ ^ r d<p r de sine |r(fi,0,0,^)|l (ii) 

471" Jo Jo Jo 

The result of numerical calculation is shown in FiglSl All response functions decrease 
above 100 MHz. This is the common feature of interferometers because the GW signal is 
destructively integrated in the cavity and is canceled out. Below 100 MHz, the response 
functions of SRI and LMI also decrease because the GW signal is partially canceled out 
during round trips of light beams in the folded cavities. On the other hand, in the case of 
FPMI, the GW signal is more constructively integrated in the cavity and is more amplified, 
as the frequency gets lower. At 100 MHz, SRI and LMI have almost the same sensitivity, 
while FPMI has sensitivity worse by a factor ^ 3.3. This is because FPMI integrates GW 
signals less efficiently than other detectors, as mentioned in Sec. Ill A[ This difference 
becomes significant when one takes the correlation of two detectors into account, which 
results in a factor of (3.3)^ ~ 10 in sensitivity to GWB energy density flg^. 
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FIG. 5: (Color online). GW response function %-msif)- Vertical lines are plotted at multiples of 
lOOMHz. 



At the end of this section, let us consider the best sensitivity of SRI experimentally 
achievable with realistic parameters, which is almost the same as that of LMI. We assume 
that the sensitivity is limited only by shot noise. The magnitude of shot noise is determined 
by laser power, the arm length and the reflectivities of mirrors, and is calculated by the 
condition that the phase shift due to GWs is equal to that of quantum noise [2J], 



(12) 



where h is the reduced Planck constant, uj is the angular frequency of laser, rj is the quantum 
efficiency of a photodetector, Jq is original laser power, and j3 = {Rp — Re)/{^ — RfRe)- 
The reason why (3 is needed in the left-hand side of the Eq.f ll2l) is that, in Eq. (1A9l) . the phase 
shift due to GWs must be converted into the amplitude of a sideband field. Substituting 
Eqs.(l5]) and ^ leads to 

1 



where we defined 



hip) 



a'{9) = \Pa{9)\ 



2h 



ReTp 



(13) 



(14) 



;i - RfRe){1 - RfRe e-^^^ 
We select L = 0.75 m so that GW signal resonates at 100 MHz, lo = 1.77 x lO^^radsec"^ 
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and r] — 1. The sensitivity achievable at 100 MHz in an ideal situation is 

with bandwidth ~ 10^ Hz, where finesse is T = 'k^/RfRe/ (1 — RfRe), which is related 
to a' by the relation T = -K^Ja'RpjTl at 100 MHz. Note that a' 1.6 x 10^ is reahzed 
with reflectivities, say, = 0.99996 and i?| = (0.99998)^. 



IV. DISCUSSIONS AND CONCLUSIONS 



In real experiments, there are many advantages and disadvantages of detectors that we 
have not considered in this paper. One of the advantages of SRI is the simplicity of a 
control system. SRI has only one degree of freedom for locking the interferometer because 
clockwise and counterclockwise lights share light paths in the cavity and the Sagnac part, 
while FPMI and LMI have three degrees of freedom in operation, which are for cavities 
in both arms and the Michelson part. Another advantage of SRI is the symmetric optical 
configuration of the cavity. This means that clockwise and counterclockwise light in the 
cavity experience the same refiectivities of mirrors. Thus, SRI is expected to have high 
tolerance to the imbalance of the reflectivities and relatively smaller laser frequency noise 
than other detectors. Considering these facts, we can conclude that SRI is the best detector. 
These issues should be investigated in more detail when one constructs real detectors. 

In this paper, we investigated the GW responses of interferometers at 100 MHz. We 
considered three designs that took advantage of the first optical resonance due to cavities 
and derived the GW response functions. As a result, SRI and LMI have almost the same 
sensitivity at 100 MHz and FPMI has sensitivity worse by a factor of 3.3. This is a significant 
difference because the sensitivity is better by a factor (3.3)^ ~ 10 in sensitivity to GWB 
energy density Jlg^ when we take the correlation between two detectors into account. When 
we select the arm length L = 0.75m, laser power /q = IW and \a'\ = 1.6 x 10^, the 
best achievable sensitivity with SRI is /i 7.8 x 10"^^ Hz""*^/^ at 100 MHz with bandwidth 
~ 2 kHz. Note that our results can also be applied to detectors at other frequency bands by 
tuning the arm length and shifting the peak of sensitivity. 
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APPENDIX A: RESPONSE FUNCTION OF SRI 

The configuration of SRI is shown in Fig{T] We call the mirror A at Xi the recycling 
mirror, the mirror B at Xi the steering mirror and the mirrors at X2 and X3 the end mirrors. 
The amplitude reflectivities and transmissivities of the steering mirror and the end mirrors 
at X2 and X3 are {Rb,Tb), (-R2,^2), (-^3,^3), respectively. Those of the recycling mirror 
are +T4) for the light incident from inside the cavity and {~Ra, +Ta) for the light 

incident from outside the cavity. The angular frequency of light is uj and the arm length is 
L. We define r = L/c. Electric fields at the recycling mirror are defined in Figj6]and are 
related by the following equations, 

C, = RaD^ + TaA,, (A1) 

Cr = RADr + TAAr , (A2) 

B, = TaD, - RaA, , (A3) 

Br = TADr - RAAr . (A4) 

We assume that there is no displacement noise, for example, thermal noise, seismic noise, ra- 







■'Cr 








< 

Br 

> 







FIG. 6: Electric fields of SRI at the recycling mirror. 
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diation pressure noise etc., at 100 MHz. Then, D field after circulating the cavity experiences 
the phase shift Aujt and the GW signal S(f>{t), 



D,{t) = R,Ci{t - At) e 
Dr{t) = RcCrit-Ar) e 



i[iujT+5(pi (t)] 



(A5) 
(A6) 



where we defined the composite reflectivity of mirrors Rc = RBR2R3 and 

5Mt) = H2i{t) + (5012 (t - r) + 5031 (t - 2r) + (5013 (t - 3r) , 
50, (t) = 5031 (t) + 5013 (t - r) + 5021 (t - 2r) + 50i2(t - 3r) . 



(AT) 
(A8) 



For example, 502i denotes the phase shift due to GW when light travels from X2 to Xi in 
FigHJ This can be calculated using Eq. ([T]). 

Equations ( ]A2I) - (1A6I) can be solved separately for right-handed and left-handed fields. 
The input-output relation of the cavity for the left-handed fields becomes 



B^it) = -RAAeit) + J2TlRcRA~^Mt - 4A;r)e^*"" exp 



k=l 



^50,(t-4(fc'-i; 



k'=l 



We assume that the cavity is in resonance in the absence of GW, that is Ai{t) = Ai{t — At). 
Then, 



B,{t) = -RAAeit) \ 1 - J^TlRcRl'^ exp 



k=l 



J2sMt-<k'-i)T) 



k'=l 



Using the approximation that the GW signal is small {S(f){t) ^ 1), we obtain 

k 



Be{t) 



Ra — Rc 
1 — RaRc 



Ae,{t) exp 



k=l 



k'=l 



Therefore, the phase shift 5$^ of left-handed light caused by GW is 

) k 



S^iit) 



1 — RaRc 



Ra — Rr 



Y.T1RcRa' $^50,(t-4(fc'-l)r) 



(A9) 



(AlO) 



k=l 



k'=l 
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Fourier transforming Eq. (lAlOp and using Eqs. (1A7P and (jl]) give 

5$, = a{n,R„RA)S^e , (AH) 



in 



1 - (e^ • ni2)^ 
1 - (e^ • ni3)^ 

(A12) 



The GW signal for the right-handed hght can be obtained by simply changing the sub- 
scripts 2^3 because of the symmetry of the system. Therefore, the output of the detector 
is 

5~^aii = S~^r - S'^e 

= a{Q,RA,Rc)S<P,a (AM) 



ni3 (X) ni3 



:{i sinnr + (e, ■ ni3)(e-*^"''^"i^ - cos fir)} 



1 - (e. ■ ni3^2 
1 - (e^ ■ ni2 



^{i sin fir + (e, ■ 1112) (e-*''™^-"^^ - cos fir)} ] . (A15) 



APPENDIX B: RESPONSE FUNCTION OF FPMI 

The configuration of FPMI is shown in FiglJl The amplitude reflectivities and trans- 
missivities of end mirrors at X^;! and are {Re,Te), and of front mirrors at Xp^i and 
Xir2 are {+Rf, +Tf) for the light incident from inside the cavities and {—Rp, +Tp) for the 
light incident from outside the cavities. The arm length is L. Electric fields at the front 
mirror are defined in Fig 13 First, we will consider only one FP cavity and calculate the 
input-output relation. At the end of our calculation, we will derive the full output of FPMI. 

The fields are related by the following equations, 

Eout = —RpEin + TfEb , (Bl) 
Ea = RfEb + TpEin . (B2) 
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Ea is reflected at the end mirror and returns to the front mirror after experiencing the 
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FIG. 7: Electric fields of FPMI. 



phase shift 2r and the modulation due to GW. The relation between Ea and Eb is 

EBit) = REEA{t - 2r) e^y^^^+Hca.m ^ (^B3) 
5(t)cav{t) = 6(j)EF{t) + S(j)pE{t - t) , (B4) 

where S(f)EF denotes the phase shift due to GW when light propagates from the end mirror 
to the front mirror and the speciflc form is given by Eq. ([1]). Similarly, 6(f)FE is the one from 
the front mirror to the end mirror. From Eqs. (]Bl|) - (]B3|) . the input-output relation for the 
FP cavity becomes 

oo k 
Eoutit) = -RpEinit) + ReTI Y^^REREf-^Einit) exp[2 ^^cav{t - 2{k' - 1)t)] 

k=l k'=l 



Re - Re / ,N 

-Ein{t) exp 



1 — ReRe 



^^ '1-ReRe) " 

- k=i k'=i 



.ReTAI - ReRe_j ^ s<PUt - 2{k' - 1 



where we used the approximation Scpcavit) 1- Therefore, the phase shift 5$ caused by 
GW is 

^^^^^ ^ _ REn{l-RFRE) ^ _ 2(k' - l)r) . (B5) 

rCp — TIE 

^ k=l k'=l 

Fourier transforming Eq. (]B5|) and using Eqs. (]B4p and (jlj) gives 



^n) = ain,RE,RE)S(fcavi^) , (B6) 



in ^ ^ 1 — fe^ ■ n) 

V 

'iQ,Te, n 



X [2 sin^]r + (e, ■n)(e-*''"'=^■"-cos^]r)] , (B7) 
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where n = (X^; — ^f)/L. This formula is consistent with the previous result, Eq.(6) in 32| . 
except for an overall constant factor. 

Using the result obtained above, we can easily obtain the full-output signal for the FPMI. 
We fix the subscripts 1 and 2 to distinguish north and east arms in Figl21 and define rij = 
(X^Ei ~ ^Fi) /L, i = 1,2. For simplicity, we assume that the two front mirrors on both arms 
are located at the same place, that is, Xi? = = Xf2- This assumption is valid because 
it hardly affects the GW signal. Then, total output of FPMI is 

= aiQ,RF,RE)S<P,a (B9) 



in 



p 

ni ® 111 



\2 



{ i sinr^r + (e^ ■ ni)(e-^^^^^-"i - cosr^r) } 



1 - (e^ ■ ni 

-T^T^^ { I sinfir + (e, ■ n2)(e-*^^''-"^ - cosier) } ] . (BIO) 
APPENDIX C: RESPONSE FUNCTION OF LMI 

The configuration of LMI is shown in Figl31 The amplitude reflectivities and transmis- 
sivities of end mirrors at X^ji and are {Re,Te), and of front mirrors at X^^i and 
are {+Re, +Te) for the light incident from inside the cavities and {—Re, +Te) for the light 
incident from outside the cavities. The two mirrors at ^ci and Xc2 are completely reflec- 
tive. The arm length is L. Electric fields at the front mirror are defined in FiglHl First, we 
will consider only one L-shaped cavity and calculate the input-output relation. At the end 
of our calculation, we will derive the full output of LMI. 

The relation between the fields is the same as Eqs. fIBip and flB2p for FPMI. The relation 
between Eb and Ea is almost the same as FPMI. However, only differences are that the 
round-trip time of LMI in the cavity is 4r and that light is reflected by the mirror at Xc 
during the trip. Thus, it is given by 

Esit) = i?sE^(t-4r) e^[^""+'''^-^W] , (CI) 
6(f),av{t) = Hcpit) + 5(t)Ec{t -r)+ 6(f)cE{t - 2r) + 6(j)Ec{t - 3r) , (C2) 

where S(f)ij, I, J = F,E,C denotes the phase shift due to GW when the light propagates 
from the I-th mirror to the J-th mirror. From Eqs. (IBip . (1B2p and (ICip . the input-output 
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Eb 
-Xp 

Eout 



BS 



FIG. 8: Electric fields of LMI. 



relation for the L-shaped cavity becomes 

oo k 
Eoutit) = -RpEinit) + ReTI Y^i^F^Ef'^Eir^it) exp[^ ^^cav{t - A{k' - 1)t)] 



k=l 



Rp — Ri 
1 — RpRp 



-Ein{t) exp 



—I 



ReTp{1 — RpRp) 
Rp — Rp 



k'=l 
oo 



k=l k'=l 

where we used the approximation Scpcavit) 1- Therefore, the phase shift 5$ caused by 
GW is 



ReTp{1 — RpRp) 
Rp — Rp 



YiRpRE)"-' Yl ^'^--(^ - - 1)^) ' (C3) 



fc=l k'=l 
Fourier transforming Eq. (1C3P and using Eqs. (lC2p and (jl]) gives 

6Mn) = a{n,Rp,RE)6<f,,,{n) , 



(C4) 



UJ 

in 



X 



ncE (K) ncE 
1 - (e^ ■ nc£;)2 

ncp <S) ncp 
1 - (e^ ■ ucpY 



{i sin fir + (e, ■ nci?)(e-*^"'=^ "c^^ - cos fir)} 
■ { i sinfir(l - 2e"^^^^^-"^^ cos fir) 



COS 2fir — cos fir) } 



a{Q,Rp,RE) = 



ReTp 



{Rp - Re){1 - RpRp e-^^-) 
where ncp = (^p - Xc)/LAncp = (Xp - Xc)/L. 



(C5) 
(C6) 
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Using the result obtained above, we can easily obtain the full-output signal for the LMI. 
We fix the subscripts 1 and 2 in order to distinguish the two arms in Figj3l and define 
111 = nc2E = —^ciF and n2 = —'nc2F = ^ciE- For simplicity, we assume that the two 
front mirrors and the two end mirrors on both arms are located at the same place, that is, 
X.p = = X.p2 and X.p = X.pi = respectively. This assumption is valid because it 
hardly affects the GW signal. Then, total output of LMI is 

6~^aii = 5~$i - S~(^2 

= a{n,Rp,RE)6\u (C7) 



p 

n2 ® n2 



{ i sin fir {e'^P^ + e'^^^ - 2 cos fir) 



1 - (e. ■ n2)2 

+ (e, ■ n2)(e"^(^^+P^) + cos2fiL - (e"*^^ + e~'P^) cos fir) } 

111 (g) 111 



{ i sin fir (e"*^^ + e-'^^ - 2 cos fir) 



1 - (e, ■ ni)2 

+ (e, • n,){e''^P'+P'^ + cos2fir - (6"'^^ + e-'P') cos fir) } ] , 
where we defined pi = fir(e2 ■ iii) and p2 = fir(e^ ■ n2). 
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